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ABSTRACT 

We have determined the visual orbit for the spectroscopic binary i Pegasi with 
interferometric visibility data obtained by the Palomar Testbed Interferometer in 
1997. i Pegasi is a double-lined binary system whose minimum masses and spectral 
typing suggests the possibility of eclipses. Our orbital and component diameter 
determinations do not favor the eclipse hypothesis: the limb-to-limb separation of the 
two components is 0.151 ± 0.069 mas at conjunction. Our conclusion that the i Peg 
system does not eclipse is supported by high-precision photometric observations. 

The physical parameters implied by our visual orbit and the spectroscopic orbit of 
Fekel and Tomkin (1983) are in good agreement with those inferred by other means. 
In particular, the orbital parallax of the system is determined to be 86.9 ±1.0 mas, 
and masses of the two components are determined to be 1.326 ± 0.016 Mq and 0.819 
± 0.009 M respectively. 

Subject headings: binaries: spectroscopic — stars: fundamental parameters — stars: 
individual (i Pegasi) — techniques: interferometric 
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1. Introduction 

i Pegasi (HR 8430, HD 210027) is a nearby, short-period (10.2 d) binary system with a F5V 
primary and a ~ G8V secondary in a circular orbit, l Peg was first discovered as a single-lined 
spectroscopic binary by Campbell (1899), and the first spectroscopic orbital elements were 
estimated by Curtis (1904). Several other single- line studies were made, notably Petrie and Phibbs 
(1949) and Abt and Levy (1976). In the context of a lithium abundance study, Herbig (1965) 
noted that lines from the i Peg secondary were visible at red wavelengths. Lithium abundances for 
both the primary ( Herbig 1965 , |Conti fc Danzinger 196€ , Duncan 1981 , Lyubimkov et al. 1991| ) 



and the secondary ( Fekel fc Tomkin 1983 , Lyubimkov et al. 1991 ) indicate the system is very 



young (~8x 10 7 yr, [Fekel fc Tomkin 1983| , 1.7 ± 0.8 x 10 8 yr, |Lyubimkov et al. 1991Q and both 



components are close to the zero-age main sequence. Both components of i Peg are also believed 



to have solar- type abundances (Lyubimkov et al. 1991) 



Following Herbig's implicit suggestion, Fekel and Tomkin (1983, hereafter FT) made 
radial velocity measurements of both l Peg components at 643 nm, and computed a definitive 
spectroscopic orbit and inferred a probable G8V spectral classification for the secondary. FT's 
orbit was noteworthy as it indicated that the minimum masses for the two components were very 
near the model values for the spectral types, suggesting a "reasonable prospect" for eclipses in the 
system (FT). Subsequent photometric monitoring by automated photometry projects in Arizona, 
at Palomar Observatory, and in Pasadena failed to show any evidence for eclipses (see §|5|). FT 
also questioned synchronous rotation of the secondary. However, Gray (1984), from somewhat 
higher resolution spectroscopic data, argued that both components are in synchronous rotation. 

Herein we report a determination of the i Peg visual orbit from near-infrared, long-baseline 
interferometric visibility measurements taken with the Palomar Testbed Interferometer. PTI is 
a 110-m K-band (2 - 2.4 /im) interferometer located at Palomar Observatory, and described in 
detail elsewhere flColavita et al. 1994 , Colavita et al. 1999a ). The minimum PTI fringe spacing 



is roughly 4 mas at the sky position of l Peg, allowing us to resolve this binary system. The 
procedures we have used to determine i Peg's visual orbit are similar to other visual orbits 



- 4 - 



determined for spectroscopic binaries using the Mark III Interferometer at Mt. Wilson (Pan et 
al. 1990| , [Armstrong et al. 1992a| , [Armstrong et al. 1992b] pan et al. 1992] , [Hummel 1993[ , pan et 
aL 1993 , Hummel et al, 1994 , Hummel et al. 1995| ), and the NPOI Interferometer at Anderson 



Mesa, AZ (Hummel et al. 1998). The analogy between i Peg and the short-period, small angular 
scale binaries studied in Hummel et al. (1995) and Hummel et al. (1998) is especially apt. 



2. Observations 

Pan attempted to determine a visual orbit for l Peg using the Mark III interferometer at 
Mt. Wilson, but the significant brightness difference in the two components at 800 nm made the 
observations difficult ( Pan 199~ ). The apparent contrast ratio in the i Peg system decreases in the 



K-band, allowing a reliable orbit determination with PTI observations. 

The observable used for these observations is the fringe contrast or visibility (squared) of an 
observed brightness distribution on the sky. Normalized in the interval [0,1], a single star exhibits 
visibility modulus given in a uniform disk model by: 

2MjrB9/X) 

ttBO/X { ' 

where J\ is the first-order Bessel function, B is the projected baseline vector magnitude at the 
star position, 9 is the apparent angular diameter of the star, and A is the center-band wavelength 
of the interferometric observation. (We consider corrections to the uniform disk model from limb 
darkening in §[|.) The expected squared visibility in a narrow pass-band for a binary star such as 
i Peg is given by: 

V nbW ~ (1 + r) 2 { ' 

where V\ and V2 are the visibility moduli for the two stars alone as given by Eq. [D, r is the apparent 
brightness ratio between the primary and companion, B is the projected baseline vector at the 
system sky position, and s is the primary-secondary angular separation vector on the plane of the 
sky (Pan et al. 199C, Hummel et al. 1995). The V 2 observables used in our 1 Peg study are both 



narrow-band V 2 from seven individual spectral channels ( Colavita et al. 1999a| ), and a synthetic 
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wide-band V 2 , given by an incoherent SNR- weighted average V 2 of the narrow-band channels in 
the PTI spectrometer flColavita 1999b ). In this model the expected wide-band V 2 observable is 



approximately given by an average of the narrow-band formula over the finite pass-band of the 
spectrometer: 

1 n 

= (3) 

i 

where the sum runs over the n = 7 channels with wavelengths \ covering the K-band (2 - 2.4 fim) 
of the PTI spectrometer in its 1997 configuration. Separate calibrations and hypothesis fits to the 
narrow-band and synthetic wide-band V 2 datasets yield statistically consistent results, with the 
synthetic wide-band data exhibiting superior fit performance. Consequently we will present only 
the results from the synthetic wide-band data. 

l Peg was observed by PTI on 24 nights between 2 July and 8 Sept 1997. In each night 
t Peg was observed in conjunction with calibration objects multiple times during the night. Each 
observation ("scan") was from 120 - 130 seconds in duration. For each scan we computed a 
mean V 2 value through methods described in Colavita (1999b). We assumed the measured rms 
in the internal scatter to be the error in V 2 . For the purposes of this analysis we have restricted 
our attention to four calibration objects, two primary calibrators within 5° of i Peg (HD 211006 
and HD 211432), and two ancillary calibrators within 15° of i Peg (HD 215510 and HD 217014 
-51 Pegasi). The suitability of 51 Peg (a known radial velocity variable) as a calibrator at 
PTI is addressed in Boden et al. (1998b). Table [l] summarizes the relevant parameters on the 
calibration objects used in this study. In particular we have estimated our calibrator diameters 
based on a model diameter on 51 Peg of 0.72 ± 0.06 mas implied by a linear diameter of 1.2 
± 0.1 Rq (adopted by [Marcy et al. 1997| ) and a parallax of 65.1 ± 0.76 mas from Hipparcos 
(|ESA 19971 , ferryman et al. 1997| ). 



The calibration of i Peg V 2 data is performed by estimating the interferometer system 
visibility iV 2 ys ) using calibration sources with model angular diameters, and then normalizing the 
raw l Peg visibility by V 2 ys to estimate the V 2 measured by an ideal interferometer at that epoch 
(Mozurkewich et al. 1991, Boden et al. 199"8a| ). We calibrated the i Peg V 2 data in two different 
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kJ-LVj OcjJciI CXtJAJll 


1_J lCXillL. VV -LVji 


Name 


Type 


Magnitude 


From l Peg 


Model 51 Peg 


HD 211006 


K2III 


5.9 V/3.4 K 


3.6° 


1.06 ± 0.05 


HD 211432 


G9III 


6.4 V/3.7 K 


3.2° 


0.70 ± 0.05 


HD 215510 


G6III 


6.3 V/3.9 K 


11° 


0.85 ± 0.06 


HD 217014 


G2.5V 


5.9 V/4.0 K 


12° 


(0.72 ± 0.06) 



Table 1: 1997 PTI l Peg Calibration Objects Considered in our Analysis. The relevant parameters 
for our four calibration objects are summarized. The apparent diameter values are determined by 
a fit to our V 2 data calibrated with respect to a model diameter for HD 217014 (51 Peg) of 0.72 ± 
0.06 mas ( Marcy et al. 199% [ESA 1997Q . 



ways: (1) with respect to the two primary calibration objects, resulting in our primary dataset 
containing 112 calibrated observations over 17 nights, and (2) an unbiased average of the primary 
and ancillary calibrators, resulting in our secondary dataset containing 151 observations over 24 
nights. The motivation for constructing these two datasets, which are clearly not independent, is 
that the determination of the orbital solution and component diameters is sensitive to calibration 
uncertainties. Comparison of the solutions derived from the two datasets allow us to quantitatively 
assess this uncertainty. 



Orbit Determination 



The estimation of the i Peg visual orbit is made by fitting a Keplerian orbit model with 
visibilities predicted by Eqs. |2| and [3| directly to the calibrated (narrow-band and synthetic 
wide-band) V 2 data on l Peg (see Armstrong et al. 1992b , Hummel 1993] , Hummel et al. 1995|) . 
The fit is non-linear in the Keplerian orbital elements, and is therefore performed by non-linear 
least-squares methods (i.e. the Marquardt-Levenberg method, Press et al. 1992] ). As such, this 
fitting procedure takes an initial estimate of the orbital elements and other parameters (e.g. 
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component angular diameters, brightness ratio), and refines the model into a new parameter set 
which best fits the data. However, the chi-squared surface has many local minima in addition 
to the global minimum corresponding to the true orbit. Because Marquardt-Levenberg strictly 
follows a downhill path in the x 2 manifold, it is necessary to thoroughly survey the space of 
possible binary parameters to distinguish between local minima and the true global minimum. In 
the case of t Peg the parameter space is significantly narrowed by the high-quality spectroscopic 
orbit and inclination constraint near 90° (FT). Furthermore, the Hipparcos distance determination 
sets the rough scale of the semi-major axis ( |ESA 1997| ). 

In addition, as the V 2 observable for the binary (Eqs. ||| and ||) is invariant under a rotation of 
180°, we cannot differentiate between an apparent primary /secondary relative orientation and its 
mirror image on the sky. In order to follow the FT convention for To at primary radial velocity 
maximum, in our analysis of i Peg we have defined To to be at a component separation extremum, 
yielding an extremum in component radial velocities for the circular orbit. We have additionally 
required our fit To to be within half a period of the projected FT determination to differentiate 
between primary radial velocity maximum and minimum. Even with our determination of To so 
defined there remains a 180° ambiguity in our determination of the longitude of the ascending 
node, Q. 

We used a preliminary orbital solution computed by Pan (1996) by separation vector 



techniques (see Pan et al. 199C for a discussion of the method), and refined it into the best-fit 
orbit shown here. We further conducted an exhaustive search of the binary parameter space that 
resulted in the same best-fit orbit, which is in fact the global minimum in the x 2 manifold. 

Figure [l] depicts the apparent relative orbit of the i Peg system. Most striking is the 
observation that the circular orbit of the system (see below) is very nearly eclipsing. From our 
primary dataset we find a best fit orbital inclination of 95.67 ± 0.21 degrees. With model angular 
diameters of 1.0 and 0.7 mas for the primary and secondary components respectively (§B), and 
an apparent semi-major axis of 10.33 ± 0.10 mas, this inclination is about 0.87° from apparent 
limb-to-limb contact. This is consistent with the lack of photometric evidence for eclipses despite 
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several photometry campaigns on the i Peg system (§H). 

Table ^ lists the complete set of V 2 measurements in the primary dataset and the prediction 
based on the best-fit orbit model for t Peg. Figure ^ shows two graphical comparisons between 
our V 2 data on i Peg and the best-fit model predictions. Figure |2|a gives four consecutive nights 
of PTI V 2 data from our primary dataset on l Peg (18 - 21 July 1997), and V 2 predictions based 
on the best-fit model for the system. Figure §b gives an additional seven consecutive nights (12 
- 18 August 1997) with the same quantities plotted. These are the two longest consecutive-night 
sequences in our data set. The model predictions are seen to be in excellent absolute and statistical 
agreement with the observed data, with a primary dataset average absolute V 2 deviation of 0.014, 
and a x 2 P er Degree of Freedom (DOF) of 0.75. 

Figure [3| gives two examples of the x 2 fit projected into orbital parameter subspaces. Figure 
||a shows a surface of x 2 /DOF projected into the subspace of orbit semi-major axis and relative 
component brightness, with all other parameters held to their best-fit values. Inset is a closeup 
of a contour plot of the x 2 /DOF surface indicating location of the best-fit parameter values, and 
contours at +1, +2, and +3 of x 2 /DOF significance. Figure ^b gives the x 2 /DOF surface in the 
subspace of orbital inclination and longitude of the ascending node. Again, the inset gives best-fit 
parameter values, and contours at +1, +2, and +3 of x 2 /DOF significance. All indications are 
that the best-fit model for the i Peg system is in excellent agreement with our V 2 data, and that 
data uniquely constrain the parameters of the visual orbit. 

Spectroscopic (from FT) and visual orbital parameters of the i Peg system are summarized 
in Table ||. We present the results for our primary and secondary datasets separately. For the 
parameters we have estimated from our interferometric data we quote a total one-sigma error in 
the parameter estimates, and the one-sigma errors in the parameter estimates from statistical 
(measurement uncertainty) and systematic error sources. In our analysis the dominant forms 
of systematic error are: (1) uncertainties in the calibrator angular diameters (Table |I|); (2) the 
uncertainty in our center-band operating wavelength (Aq ~ 2.2 /rai), which we have taken to be 
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Fig. 1. — Visual Orbit of i Pegasi. The relative visual orbit of l Peg is depicted, with the primary 
and secondary rendered at To (maximum primary radial velocity) and apparent conjunction. The 
inset shows a closeup of the system at apparent conjunction. By our model the i Peg orbit is nearly, 
but not quite eclipsing, being approximately 0.87° in inclination from apparent grazing eclipses. 



-10- 



-a 

o 



a! 

> 



83 
o 



> 



07/18 
00:00 
1 



Calendar Date 

07/18 07/19 07/19 07/20 07/20 07/21 07/21 07/22 
12:00 00:00 12:00 00:00 12:00 00:00 12:00 00:00 



0.8 
0.6 
0.4 
0.2 

0.1 
0.05 


-0.05 
-0.1 



a) 



i Peg Model Predict 

V 2 Data -*- 




50647 50647.5 50648 50648.5 50649 50649.5 50650 50650.5 50651 

Modified Julian Date (days) 



08/12 08/13 
00:00 00:00 
1 



Calendar Date 
08/14 08/15 08/16 
00:00 00:00 00:00 



08/17 08/18 08/19 

00:00 00:00 00:00 



0.8 
0.6 
0.4 
0.2 

0.1 
0.05 


-0.05 
-0.1 



b) 



i 


h 


& 


* 




f 








f 


I ! 


i 


sjs 


\ k 






± 

as.. 





i Peg Model Predict 

V 2 Data ^ * 




50672 50673 50674 50675 50676 50677 50678 50679 

Modified Julian Date (days) 



Fig. 2. — V 2 Fit of i Pegasi. a) Four consecutive nights (18 - 21 July 1997) of calibrated V 2 data 
on l Peg, and V 2 predictions from the best-fit model for the system. In the lower frame we give 
V 2 residuals between the calibrated data and best-fit model, b) An additional seven consecutive 
nights (12 - 18 August 1997) of data on l Peg, with model predicts and fit residuals. The model 
is in good agreement with the calibrated data, with a x 2 /DOF of 0.75 and an average absolute V 2 
residual of 0.014. 
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Fig. 3. — x 2 /DOF Fit Surfaces for t Pegasi Primary Dataset. a) x 2 /DOF surface in the subspace 
of orbit semi-major axis and relative component brightness. Inset is a closeup of a contour plot 
surface indicating location of the best-fit parameter values, and contours at +1, +2, and +3 of 
X 2 /DOF significance, b) \ 2 /DOF surface in the subspace of orbital inclination and longitude of 
the ascending node, with inset giving surface contour closeup. 
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20 nm (~1%); (3) the geometrical uncertainty in our interferometric baseline ( < 0.01%); and 
(4) uncertainties in orbital parameters we have constrained in our fitting procedure (e.g. period, 
eccentricity). Different parameters are affected differently by these error sources; our estimated 
uncertainty in the i Peg orbital inclination is dominated by measurement uncertainty, while the 
uncertainty in the angular semi-major axis is dominated by uncertainty in the wavelength scale. 
Conversely, we have assumed that all the uncertainty quoted by FT in the t Peg spectroscopic 
parameters is statistical. Finally, we have listed the level of statistical agreement in the visual 
orbit parameters in our two solutions (the absolute residual between the two estimates divided by 
the RSS of their statistical errors). The two solutions are in good statistical agreement, giving us 
confidence we have properly characterized our calibration uncertainties. 

Particularly remarkable is the agreement between To (quoted as the epoch of maximum 
primary radial velocity for the t Peg circular orbit) and period as determined by FT, and To as 
determined in our primary dataset, separated from the FT determination by 523 cycles. FT quote 
an l Peg period accurate to roughly 1 part in 10 6 , resulting in a propagated uncertainty in To at 
the epoch of our observations of 7 x 10~ 3 days. This FT-extrapolated To differs from our 1997 
To determination by 8 x 10 -4 days, an agreement of roughly 0.1 sigma. A similar comparison 
with the secondary dataset solution is less spectacular, an agreement at 0.7 sigma. Clearly the 
extraordinary quoted accuracy of the t Peg period determination by FT (made by combining their 
1977 - 1982 data with spectroscopy from the mid-30s - Petrie fc Phibbs 1949| ) seems well justified 



compared to our visual orbit. Consequently we have assumed the FT value for the i Peg period. 

Following FT we have assumed a circular orbit for the system. Fitting our primary dataset 
for an eccentricity in the system yields an estimate of 1.5 x 10~ 3 ± 1.3 x 10~ 3 . The assumption 
of a circular orbit seems well justified. 
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10.32 ± 0.11 (0.02/0.11) 


0.35 


A K (mag) 




1.610 ± 0.021 
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Table 3: Orbital Parameters for i Peg. Summarized here are the apparent orbital parameters for the 
i Peg system as determined by FT, and our PTI primary and secondary datasets. For parameters 
estimated from our PTI observations we separately quote one sigma errors from both statistical and 
systematic sources (listed as <T s tat/&sys), and the total error as the sum of the two in quadrature. 
We have also included the level of statistical agreement between visual orbit parameters from our 
two solutions; the parameters estimated separately from the primary and secondary datasets are in 
good agreement in relation to the statistical component of their error estimates. We have quoted 
the longitude of the ascending node parameter (fi) as the angle between local East and the orbital 
line of nodes (and the relative position of the secondary at To), measured positive in the direction 
of local North. Due to the degeneracy in our V 2 observable there is a 180° ambiguity in Finally, 
the fit x 2 / D F and mean absolute V 2 residual (|-Ry2|) is listed for both solutions. 
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4. Physical Parameters 

Physical parameters derived from the i Peg primary dataset visual orbit and the FT 
spectroscopic orbit are summarized in Table ||. We use the primary dataset solution because 
it is the most free from possible sky position-dependent systematic effects (as the secondary 
dataset includes the ancillary calibrators), but we note the two orbital solutions yield statistically 
consistent results. Notable among the physical parameters for the system is the high-precision 
determination of the component masses for the system, a virtue of the precision of the FT radial 
velocities on both components and the high inclination of the orbit. We estimate the masses of 
the F5V primary and putative G8V secondary components as 1.326 ± 0.016 Mq and 0.819 ± 
0.009 Mq respectively. Our mass values agree well with mass estimates of 1.33 ± 0.08 Mq and 
0.9 ± 0.2 Mq respectively made by Lyubimkov et al. (1991) based on evolutionary models and 
spectroscopic measurements of component effective temperatures and surface gravities. 

The Hipparcos catalog lists the parallax of i Peg as 85.06 ± 0.71 mas ( [ESA 1997| ) . The 
distance determination to i Peg based on the FT radial velocities and our apparent semi-major 
axis and inclination is 11.51 ± 0.13 pc, corresponding to an orbital parallax of 86.91 ± 1.0 mas, 
consistent with the Hipparcos result at roughly 2% and 1.5 sigma. 

FT list main-sequence model linear diameters for the two c Peg components as 1.3 and 0.9 Rq 
respectively (FT). At a distance of approximately 11.5 pc this corresponds to apparent angular 
diameters of 1.0 and 0.7 mas for the primary and secondary components respectively. We have 
fit for the uniform-disk angular diameter for both components as a part of the orbit estimation, 
and find best fit apparent diameters of 0.98 =L 0.05 and 0.70 ± 0.10 mas. Because we have limited 
spatial frequency coverage in our data, following Mozurkewich et al. (1991) and Quirrenbach et 
al. (1996) we have estimated the limb-darkened diameters of the components from a correction 
to the uniform-disk diameter based on the solar limb-darkening at 2 ^m given by Allen (1982). 
The limb-darkened diameters for the primary and secondary components are 1.0 ± 0.05 and 0.71 
± 0.10 mas respectively. For both the primary and secondary components our fits for apparent 
diameter are in good agreement with main-sequence model diameters. 
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The observed K-magnitude of the i Peg system (2.623 ± 0.016 - Carrasco et al. 1991| , 2.656 db 
0.002 - Bouchet et al. 1991] ) and our estimates of the distance and relative K-photometry (Table |3|) 
of the system allows the determination of the absolute magnitude of both components separately. 
Using the Bouchet et al. (1991) K-photometry we obtain M K values of 2.574 ± 0.025 and 4.182 
± 0.030 for the primary and secondary components respectively. Both of these M# values are 
consistent (within quoted scatter) to the empirical mass-luminosity relation for nearby low-mass, 
main-sequence stars given by Henry & McCarthy (1992, 1993). In particular, our value for 
the primary is 0.010 mag brighter than the mass-luminosity prediction flHenry &: McCarthy 1992 ), 
while the 4.18 Mx value for the secondary is roughly 0.28 magnitudes dimmer than the prediction 



(Henry & McCarthy 1993). Both values are well within the quoted scatter of the mass- luminosity 
models. A second check on the absolute K-magnitude estimates can be extracted from the model 
calculations of Bertelli et al. (1994), who predict absolute K-magnitudes of 2.616 ± 0.048 and 
4.254 ± 0.039 for our estimated primary and secondary masses respectively for main-sequence 
stars with solar-type abundances at an age of 1.7 ± 0.8 x 10 8 yr flLyubimkov et al. 1991|). 



5. Eclipse Search 

A critical test of our visual orbit model is a high-precision photometric search for eclipses 
in i Peg. Combined with our visual orbit (Table ||) , our measured diameters (Table |]) imply an 
apparent limb-to-limb separation at conjunction of 0.151 db 0.069 mas (using our limb-darkened 
diameter estimates). Our visual orbit and fit diameters do not favor the FT conjecture of possible 
eclipses in the i Peg system. Conversely, were the inclination of the orbit near 90°, there would be 
significant primary eclipses with a duration of a few hours (6.8 hr for i = 90° - FT), and as large 
as 0.6 mag in V-band. 

Several individuals have searched for signs of eclipses in the t Peg system. In 1997 both Van 
Buren with the 60" telescope at Palomar (1997) and one of us (C.D.K.) at the Robinson Rooftop 
Observatory at Caltech in Pasadena (Koresko 1997) searched for eclipses during primary and 



secondary eclipse opportunities respectively. Both searches resulted in non-detections at about the 
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Physical 


Primary 


Secondary 


Parameter 


Component 


Component 


a (1(T 2 AU) 


4.54 ± 0.03 (0.03/0.0002) 


7.35 ± 0.03 (0.03/0.0003) 


Mass (M ) 


1.326 ± 0.016 (0.016/0.0001) 


0.819 ± 0.009 (0.009/0.0001) 


Sp Type (FT) 


F5V 


G8V 


Model Diameter (mas) 


1.0 


0.7 


UD Fit Diameter (mas) 


0.98 ± 0.05 (0.01/0.05) 


0.70 ± 0.10 (0.03/0.10) 


LD Fit Diameter (mas) 


1.0 ± 0.05 (0.01/0.05) 


0.71 ± 0.10 (0.03/0.10) 


System Distance (pc) 


11.51 ± 0.13 (0.05/0.12) 


ir orb (mas) 


86.91 ± 1.0 (0.34/0.94) 


M K (mag) 


2.574 ± 0.025 (0.010/0.024) 


4.182 ± 0.030 (0.019/0.028) 



Table 4: Physical Parameters for t Peg. Summarized here are the physical parameters for the i Peg 
system as derived from the orbital parameters in Table |3|. As for our PTI-derived orbital parameters 
we have quoted both total error and separate contributions from statistical and systematic sources 
(given as <J sm /a sys ). 
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0.1 mag levels. 

More comprehensive and sensitive than the Southern California searches has been the 
program conducted by the Automated Astronomy Group at Tennessee State University, i Peg was 
observed photometrically in 1984 with the Phoenix-10 automatic photoelectric telescope (APT) in 
Phoenix, AZ, and again in 1997-98 with the Vanderbilt/Tennessee State 16-inch APT at Fairborn 
Observatory near Washington Camp, AZ, in order to search for possible eclipses suggested by FT. 
Both telescopes observed i Peg once per night through a Johnson V filter with respect to the 
comparison star HR 8441 (HD 210210, Fl IV) in the sequence C,V,C,V,C,V,C, where C is the 
comparison star and V is i Peg. Three differential magnitudes (in the sense V-C) were computed 
from each nightly sequence, corrected for differential extinction, and transformed to the Johnson 
system. The three differential magnitudes from each sequence were then averaged together and 
treated as single observations thereafter. Because of the lack of accurate standardization in the 
Phoenix-10 data set, a -0.027 mag correction was added to each observation to bring those data 
in line with the 16-inch observations. The observations are summarized in Table |5[ Column 4 
gives the standard deviation of a single nightly observation from the mean of the entire data set 
and represents a measure of the precision of the observations. Further details on the telescopes, 
data acquisition, reductions, and quality control can be found in Young et al. (1991) and Henry 
(1995a,b). 

The photometric observations summarized in Table |H| are plotted in Figure |I] against orbital 
phase of the binary computed from the FT-defined Tq and period. For inclinations allowing 



APT 


JD Range 


# Obs. 


Std. Dev. 




(+2400000) 




(mag) 


10- inch 


45703 - 46065 


78 


0.0109 


16-inch 


50718 - 50829 


66 


0.0032 



Table 5: Summary of APT Photometry on i Peg. 
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eclipses of the two components, the phases of conjunction coinciding with primary and secondary 
eclipse opportunities are 0.25 and 0.75 respectively. FT estimated the total duration of a central 
eclipse (i = 90°) to be roughly 6.8 hours or 0.027 phase units. Our photometric observations 
exclude this possibility and show no evidence for any partial eclipse to a precision of around 0.003 
mag. The time of conjunction is uncertain by no more than a few minutes, and gaps in the data 
around the time of conjunction are no larger than about 0.005 phase units (1.2 hours). Thus, 
the possibility of all but the briefest of grazing eclipses are excluded by the APT photometry. In 
particular, using the two points nearest the primary conjunction opportunity (at -1.29 and +1.22 
hours relative to the predicted conjunction respectively) constrain |90 — i\ to be greater than 4.07° 
and 4.10° respectively at greater than 99% confidence, based on the model diameters and M v 
estimates of 3.4 and 5.8 for the primary and secondary components respectively. 

The components of most close binaries with orbital periods less than about one month rotate 



synchronously with the orbital period due to tidal action between the components (e.g. Fekel & 



Eitter 1989). Such synchronous rotation is expected in l Peg and is confirmed by the rotational 



broadening measurements of FT and Gray (1984) (c.f. Wolff & Simon 1997). If the G8V secondary. 



which is much more convective than the F5V primary, is rotating synchronously, it would be 
expected to be photometrically variable on the orbital period at the level of a few percent due 
to starspot activity ( |Henry et al. 1999 ). In fact, i Peg is listed as a suspected variable star by 



Petit (1990), who reports variability at the 0.02 mag level in V. FT estimate that the secondary is 
roughly 2.7 mag fainter in the V band than the primary, so any apparent photometric variability 
of the secondary component will be diluted by a factor of about 12 by the primary component. 

In order to search for this possible photometric variability in i Peg, we performed a 
periodogram analysis of the 16-inch APT data. The analysis reveals a photometric period that is 
identical, within its uncertainty, to the spectroscopic period, a result that is consistent with the 
assumption of synchronous rotation. Likewise, the amplitude of 0.0037 mag, scaled by a factor 
of 12, results in a 4.4% variation, similar to the variability expected from rotational modulation 
of the spotted surface of the secondary diluted by the emission of the primary. Based on these 
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Conj Offset (hr) 
-6 -3 3 6 



-2.36 
-2.34 
-2.32 
-2.3 
-2.28 
-2.26 



84 Phoenix- 10 APT 
97/98 16-inch APT 




A- 



T T >A • -f , 

A A # A A A 



Primary Conjunction 



-2.295 



0.22 0.25 

S'iT A | A A A 

A T * T , A Vi* T . _ 



Secondary Conjunction 



0.28 



0.25 0.5 0.75 

I Peg Phase (dimensionless) 



Fig. 4. — Photometric Observations of i Peg. Differential photometric observations of i Peg from 
the Phoenix-10 APT (open triangles) and the Vanderbilt /Tennessee State University 16-inch APT 
(filled triangles) plotted against orbital phase of the binary computed following FT. Phase 0.25 
represents a time of conjunction with the secondary in front (primary eclipse opportunity). Inset 
we show a closeup of the data around the primary eclipse opportunity. (We have added a second 
horizontal scale relative to the eclipse opportunity in units of hours; a full eclipse in the i Peg system 
would be roughly 7 hours in duration.) The photometric observations exclude the possibility of all 
but the briefest of grazing eclipses in the l Peg system. 
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results, we conclude that i Peg is a low-amplitude variable star. 

6. Summary 

We have presented the visual orbit for the double- lined binary system i Pegasi, and derived 
the physical parameters of the system by combining it with the earlier spectroscopic orbit of Fekel 
and Tomkin. The derived physical parameters of the two young stars in i Peg are in reasonable 
agreement with the results of other studies of the system, and theoretical expectations for stars 
of these types. Noted by FT, the l Peg system is nearly eclipsing; because our model visual orbit 
is so close to producing observable eclipses we have further presented high-precision photometric 
data which is consistent with our visual orbit model. 

l Peg represents a prototype of the binary system that PTI is well-suited to measure; the 
large magnitude difference between components in the visible is significantly mitigated in the 
near-infrared, making the accurate determination of the system parameters feasible. 

Part of the work described in this paper was performed at the Jet Propulsion Laboratory, 
California Institute of Technology under contract with the National Aeronautics and Space 
Administration. Interferometer data was obtained at the Palomar Observatory using the NASA 
Palomar Testbed Interferometer, supported by NASA contracts to the Jet Propulsion Laboratory. 

Automated astronomy at TSU has been supported for several years by the National 
Aeronautics and Space Administration and by the National Science Foundation, most recently 
through NASA grants NCC2-977 and NCC5-228 (which supports TSU's Center for Automated 
Space Science) and NSF grants HRD-9550561 and HRD-9706268 (which supports TSU's Center 
for Systems Science Research). 

We wish to thank the anonymous referee for his many positive contributions to the accuracy 
and quality of this manuscript, and his forbearance in the review process. 
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TABLE 2 

i Peg Primary Dataset, Best-Fit Model Predictions, and Residuals 



Mod Julian Date 


Calendar Date/UT 


Measured V 2 


Error 


V 2 Model Pred 


V 2 Resid 


50647.30181 


7/18/97 7:14:36 


0.846 


0.038 


0.815 


0.031 


50647.31556 


7/18/97 7:34:23 


0.819 


0.043 


0.835 


-0.016 


50647.32495 


7/18/97 7:47:55 


0.872 


0.030 


0.846 


0.026 


50647.33541 


7/18/97 8 


02:59 


0.849 


0.036 


0.855 


-0.007 


50647.34728 


7/18/97 8 


20:04 


0.866 


0.022 


0.862 


0.004 


50647.35482 


7/18/97 8 


30:56 


0.880 


0.047 


0.865 


0.015 


50647.36287 


7/18/97 8 


42:31 


0.895 


0.026 


0.866 


0.029 


50647.37229 


7/18/97 8 


56:06 


0.865 


0.036 


0.866 


-0.001 


50647.37981 


7/18/97 9:06:55 


0.865 


0.024 


0.865 


0.000 


50647.38748 


7/18/97 9:17:58 


0.830 


0.030 


0.863 


-0.032 


50647.39693 


7/18/97 9:31:34 


0.841 


0.021 


0.859 


-0.018 


50647.40395 


7/18/97 9:41:41 


0.841 


0.015 


0.855 


-0.014 


50647.4113 


7/18/97 9:52:15 


0.830 


0.028 


0.851 


-0.021 


50647.42058 


7/18/97 10:05:38 


0.828 


0.019 


0.846 


-0.018 


50647.42789 


7/18/97 10:16:09 


0.832 


0.020 


0.841 


-0.009 


50648.32444 


7/19/97 7 


47:11 


0.722 


0.035 


0.691 


0.031 


50648.33364 


7/19/97 8 


00:26 


0.634 


0.025 


0.677 


-0.043 


50648.34342 


7/19/97 8 


14:31 


0.664 


0.028 


0.663 


0.001 


50648.35301 


7/19/97 8 


28:20 


0.662 


0.028 


0.648 


0.014 


50648.36952 


7/19/97 8 


52:06 


0.657 


0.026 


0.624 


0.032 


50648.37877 


7/19/97 9 


05:25 


0.599 


0.016 


0.612 


-0.013 


50648.38843 


7/19/97 9 


19:20 


0.613 


0.021 


0.598 


0.015 


50648.39773 


7/19/97 9 


32:43 


0.585 


0.019 


0.586 


-0.000 


50649.36696 


7/20/97 8 


48:25 


0.409 


0.019 


0.405 


0.004 


50649.38384 


7/20/97 9 


12:43 


0.368 


0.015 


0.369 


-0.002 


50649.30416 


7/20/97 9 


41:59 


0.356 


0.014 


0.353 


0.003 


50649.4217 


7/20/97 10:07:14 


0.344 


0.013 


0.359 


-0.015 


50649.43655 


7/20/97 10:28:37 


0.369 


0.017 


0.375 


-0.006 


50650.34729 


7/21/97 8 


20:05 


0.851 


0.018 


0.842 


0.009 


50650.37801 


7/21/97 9:04:19 


0.790 


0.015 


0.800 


-0.010 


50650.40195 


7/21/97 9:38:48 


0.691 


0.012 


0.710 


-0.019 


50650.41566 


7/21/97 9:58:32 


0.673 


0.075 


0.657 


0.016 


50650.41948 


7/21/97 10:04:02 


0.592 


0.031 


0.646 


-0.055 


50659.33293 


7/30/97 7:59:25 


0.774 


0.018 


0.775 


-0.001 


50659.34011 


7/30/97 8 


09:45 


0.805 


0.022 


0.753 


0.052 


50659.36826 


7/30/97 8 


50:17 


0.656 


0.024 


0.657 


-0.001 


50659.37567 


7/30/97 9:00:57 


0.649 


0.017 


0.632 


0.017 


50659.40534 


7/30/97 9:43:41 


0.554 


0.017 


0.543 


0.011 


0UD0y.41zZ0 


7/30/97 9:53:38 


u.04y 


n ni7 


U.ozo 


n no/i 

U.UZ4 


50659.44485 


7/30/97 10:40:35 


0.479 


0.018 


0.457 


0.022 


50659.45221 


7/30/97 10:51:10 


0.453 


0.015 


0.445 


0.008 


50659.48217 


7/30/97 11:34:19 


0.434 


0.024 


0.411 


0.023 


50659.48931 


7/30/97 11:44:36 


0.427 


0.025 


0.407 


0.020 


50660.31914 


7/31/97 7 


39:33 


0.786 


0.022 


0.783 


0.003 


50660.3346 


7/31/97 8 


01:49 


0.850 


0.018 


0.829 


0.021 


50660.35318 


7/31/97 8 


28:34 


0.839 


0.029 


0.845 


-0.006 


50660.36865 


7/31/97 8 


50:51 


0.857 


0.028 


0.831 


0.026 


50660.39155 


7/31/97 9 


23:49 


0.779 


0.033 


0.775 


0.004 
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TABLE 2— Continued 



Mod Julian Date 


Calendar Date/UT 


Measured V 2 


Error 


V 2 Model Pred 


V 2 Resid 


50660.40886 


7/31/97 9:48:45 


0.716 


0.028 


0.722 


-0.006 


50660.43543 


7/31/97 10:27:01 


0.619 


0.039 


0.643 


-0.023 


50660.45022 


7/31/97 10:48:18 


0.606 


0.019 


0.609 


-0.003 


50660.46569 


7/31/97 11:10:35 


0.612 


0.031 


0.585 


0.027 


50660.48055 


7/31/97 11:31:59 


0.577 


0.022 


0.565 


0.011 


50660.49702 


7/31/97 11:55:42 


0.533 


0.048 


0.555 


-0.021 


50660.50991 


7/31/97 12:14:16 


0.474 


0.033 


0.549 


-0.075 


50661.31235 


8/1/97 7:29:47 


0.871 


0.027 


0.837 


0.033 


50661.32807 


8/1/97 7:52:25 


0.850 


0.030 


0.842 


0.008 


50661.34898 


8/1/97 8:22:31 


0.830 


0.026 


0.805 


0.024 


50661.36423 


8/1/97 8:44:29 


0.784 


0.019 


0.756 


0.028 


50661.3849 


8/1/97 9:14:15 


0.694 


0.014 


0.678 


0.016 


50661.39919 


8/1/97 9:34:50 


0.638 


0.015 


0.626 


0.012 


50661.4217 


8/1/97 10:07:15 


0.570 


0.013 


0.558 


0.012 


50661.43594 


8/1/97 10:27:44 


0.538 


0.011 


0.526 


0.012 


50661.4566 


8/1/97 10:57:30 


0.498 


0.014 


0.495 


0.003 


50661.47175 


8/1/97 11:19:19 


0.478 


0.022 


0.481 


-0.003 


50661.498 


8/1/97 11:57:07 


0.472 


0.018 


0.473 


-0.002 


50672.34939 


8/12/97 8:23:07 


0.542 


0.012 


0.544 


-0.002 


50672.36492 


8/12/97 8:45:28 


0.586 


0.012 


0.588 


-0.001 


50672.41048 


8/12/97 9:51:05 


0.695 


0.030 


0.682 


0.013 


50672.42547 


8/12/97 10:12:40 


0.702 


0.025 


0.699 


0.003 


50672.47332 


8/12/97 11:21:35 


0.699 


0.017 


0.704 


-0.005 


50672.48845 


8/12/97 11:43:22 


0.696 


0.018 


0.692 


0.004 


50673.29096 


8/13/97 6:58:59 


0.336 


0.010 


0.340 


-0.004 


50673.33522 


8/13/97 8:02:42 


0.344 


0.010 


0.349 


-0.006 


50673.37336 


8/13/97 8:57:37 


0.360 


0.007 


0.357 


0.003 


50673.40543 


8/13/97 9:43:49 


0.352 


0.017 


0.358 


-0.007 


50674.30151 


8/14/97 7:14:10 


0.594 


0.011 


0.589 


0.005 


50674.33068 


8/14/97 7:56:10 


0.684 


0.009 


0.668 


0.016 


50674.34999 


8/14/97 8:23:58 


0.722 


0.010 


0.714 


0.007 


50674.38467 


8/14/97 9:13:55 


0.798 


0.013 


0.780 


0.018 


50675.27741 


8/15/97 6:39:28 


0.488 


0.007 


0.488 


-0.001 


50675.29329 


8/15/97 7:02:20 


0.567 


0.011 


0.569 


-0.002 


50675.30871 


8/15/97 7:24:32 


0.640 


0.011 


0.648 


-0.008 


50675.32331 


8/15/97 7:45:34 


0.706 


0.012 


0.715 


-0.008 


50675.34494 


8/15/97 8:16:42 


0.782 


0.013 


0.788 


-0.005 


50676.27174 


8/16/97 6:31:18 


0.801 


0.054 


0.841 


-0.040 


50676.3422 


8/16/97 8:12:45 


0.565 


0.011 


0.580 


-0.016 


50676.36442 


8/16/97 8:44:45 


0.496 


0.012 


0.502 


-0.006 


50676.37941 


8/16/97 9:06:21 


0.452 


0.016 


0.462 


-0.011 


50676.39791 


8/16/97 9:32:59 


0.413 


0.010 


0.428 


-0.014 


50677.2532 


8/17/97 6:04:36 


0.404 


0.013 


0.395 


0.009 


50677.29427 


8/17/97 7:03:45 


0.552 


0.020 


0.546 


0.006 


50677.33911 


8/17/97 8:08:19 


0.700 


0.018 


0.713 


-0.013 


50677.35929 


8/17/97 8:37:22 


0.760 


0.015 


0.764 


-0.004 


50678.2794 


8/18/97 6:42:19 


0.500 


0.021 


0.482 


0.018 


50678.39388 


8/18/97 9:27:10 


0.617 


0.066 


0.575 


0.042 
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TABLE 2— Continued 



Mod Julian Date 


Calendar Date/UT 


Measured V 2 


Error 


V 2 Model Pred 


V 2 Resid 


50678.46302 


8/18/97 11:06:44 


0.609 


0.075 


0.545 


0.065 


50681.28315 


8/21/97 6:47:44 


0.792 


0.023 


0.787 


0.005 


50681.32948 


8/21/97 7:54:27 


0.593 


0.014 


0.587 


0.006 


50681.35704 


8/21/97 8:34:07 


0.483 


0.013 


0.491 


-0.007 


50684.2004 


8/24/97 4:48:34 


0.437 


0.022 


0.452 


-0.015 


50684.21621 


8/24/97 5:11:20 


0.444 


0.019 


0.423 


0.021 


50684.23825 


8/24/97 5:43:04 


0.395 


0.020 


0.390 


0.005 


50684.27121 


8/24/97 6:30:32 


0.329 


0.012 


0.356 


-0.027 


50684.31539 


8/24/97 7:34:09 


0.346 


0.011 


0.338 


0.009 


50684.33993 


8/24/97 8:09:30 


0.335 


0.009 


0.339 


-0.004 


50684.35529 


8/24/97 8:31:36 


0.341 


0.007 


0.343 


-0.002 


50684.39971 


8/24/97 9:35:34 


0.363 


0.010 


0.361 


0.002 


50684.43043 


8/24/97 10:19:49 


0.381 


0.014 


0.375 


0.005 


50699.30539 


9/8/97 7:19:45 


0.890 


0.023 


0.879 


0.011 


50699.35183 


9/8/97 8:26:37 


0.843 


0.020 


0.878 


-0.035 


50699.38659 


9/8/97 9:16:41 


0.887 


0.023 


0.871 


0.016 
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